Propagating waves along spicules 
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ABSTRACT 



Alfvenic waves are thought to play an important role in coronal heating 
and acceleration of solar wind. Here we investigated the statistical properties 
of Alfvenic waves along spicules (jets that protrude into the corona) in a po- 
lar coronal hole using high cadence observations of the Solar Optical Telescope 
(SOT) onboard Hinode. We developed a technique for the automated detection of 
spicules and high-frequency waves. We detected 89 spicules, and found: (1) a mix 
of upward propagating, downward propagating, as well as standing waves (occur- 
rence rates of 59%, 21%, and 20%, respectively). (2) The phase speed gradually 
increases with height. (3) Upward waves dominant at lower altitudes, standing 
waves at higher altitudes. (4) Standing waves dominant in the early and late 
phases of each spicule, while upward waves were dominant in the middle phase. 

(5) In some spicules, we find waves propagating upward (from the bottom) and 
downward (from the top) to form a standing wave in the middle of the spicule. 

(6) The medians of the amplitude, period, and velocity amplitude were 55 km, 45 
s, and 7.4 km s _1 , respectively. We speculate that upward propagating waves are 
produced near the solar surface (below the spicule) and downward propagating 
waves are caused by reflection of (initially) upward propagating waves off the 
transition region at the spicule top. The mix of upward and downward prop- 
agating waves implies that exploiting these waves to perform seismology of the 
spicular environment requires careful analysis and may be problematic. 

Subject headings: waves — Sun: chromosphere — Sun: transition region 
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Introduction 



Coronal heating and the acceleration of the solar wind are unsolved problems in solar 
physics. The propagation of Alfven waves along magnetic field lines is one of the candidate 
mechanisms to carry energy to large distances from the surface and heat the coronal plasma. 



excited in the photosphere or chromosphere (e.g., 


Hollwcg et al 


1982; 


Tu & Marsch 


1997 


Axford et al. 


1999; 


Suzuki & Inutsuka 


2005; 


Cranmer et al. 


2007 


Antolin & Shibatal 


2010) 



It has proven difficult to show from observations that waves are driving coronal heat- 
ing. It is clear that three components should be found: the generation, propagation, and 
dissipation of waves. Low-frequency transverse waves propagating in the cor ona have been 



reported using line-qf-sight Dopplergrams from groun d-based coronagraphs (ISakurai et al. 



2002 



2010 



Tomczvk et al.l2 007 



Mcintosh et al 



2011 



Tomczvk fc Mclntoshll2009l ) and spectrometers in orbit (IBanerjee et al. 



and references therein). Here we focus on waves along chromo- 
spheric features embedded in the corona because chromospheric observations can be obtained 
at high enough spatial resolution to directly resolve the transverse motion of field lines as 
a result of Alfven waves. While there have been earlier reports (fr om ground-based data) 
of th e presence of transverse waves along spicules (for a review, see IZaaarashvili fc Erdelyi 
20091 ). filtergraph observations with the Hinode satellite (IKosugi et al.ll2007l) for the first time 
clearly resolved small sc ale oscillations in prominences ( jokamoto et al. 2007 ) and spicules 
( IDe Pontieu et al.ll2007al ). These oscillations were identified as low- frequency Alfvenic waves 
in the coronal volume. In the latter cases, it is not clear whether these waves are actually 
propagating waves or standing waves, because t he phase d i fferen ce of the oscillations along 
the field lines were not resolved. More recently, iHe et al.l ( 120091 ) introduced 4 examples of 



spicules that had high-frequency upward propagating waves with very small amplitudes. 

The purpose of this paper is to determine the properties of transverse waves for a 
large statistical sample of spicul es. We obtained observations with high ca d ence, using the 



Solar Optical Telescope ( SOT; iTsuneta et al.l 120081 ; ISuematsu et al.l 120081 ; llchimoto et al. 



2008 



Shimizu et al.l 120081 ) onboard Hinode, to resolve the phase difference of waves along 
spicules. We developed a technique for the automatic detection of spicules and investigated 
the observed waves. We present an analysis of the waves, and discuss their properties. 



2. Observation and Data Reduction 



The Hinode satellite observed a coronal hole (CH) boundary at the south pole with a 
cadence of 1.6 s from 11:06 UT to 11:59 UT on January 29, 2011. We obtained 1995 images 
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of off-limb spicules with the Ca II H-line filter (3968A, 3A bandwidth) of the SOT. The field 
of view is 55"x55"(512x512 pixel 2 ). 

Figure [TJ shows a Ca II H-line image at 11:06 UT and simultaneous coronal images 
taken with the 171A filter of the Atmospheric Imaging Assembly (AIA) onboard the Solar 
Dynamics Observatory. While it is difficult to determine the boundaries of the CH because 
of line-of-sight confusion, the connectivity of the field lines visible off-limb (revealed by deep 
exposures and image processing of 171 A images) indicates that the observed spicules are 
most likely located at the footpoints of open field lines. 

To automatically detect a spicule and track it in time, we developed a technique based on 
the following steps: (1) A radial density filter is applied to enhance visibility of faint features. 
(2) Each snapshot is averaged with 9 cons ecutive images to increase the signal-to-noise. (3) 



A high-pass filter (the a trous algorithm: IStarck et al.lll998l ) is applied to emphasize bright 



features. (4) The IDL procedure "thin. pro" is used to determine the coordinates of the 
medial axis of the bright features. (5) To qualify as a spicule, each bright feature is required 
to be more than 4" (40 pixels) long and have a deviation from a straight line that is less 
than 5° over every 4". (6) The bright feature has to contain at least 30 pixels for which the 
horizontal displacement between consecutive snapshots (1.6 s apart) is within 1 pixel (at the 
same height). If there is no spicule detected in the timestep following the current timestep, 
the condition is applied to the medial axis of the next timestep that has a defined spicule 
(as long as it is within three timesteps of the current timestep). (7) Each spicule is required 
to be longer than 8" at its maximum extent, and live longer than 40 s (25 timesteps). (8) 
Each coordinate of the medial axes is calculated at the sub-pixel level by fitting a Gaussian. 
If the center axis deduced from the fitting is 4 pixels removed from the original medial axis, 
the Gaussian fit is not used. 

Using these procedures, we have obtained timeseries' of 96 different spicules. We visually 
inspected all spicules and rejected 7 that have irregular shapes, typically because multiple 
spicules are superposed along the line of sight. After removing these irregular spicules, we 
analyze the remaining 89 spicules (online animation A). While the detection method does 
not automatically classify the spicules as type I or II, most of the observed spicules are likely 



so-called type-II spicules ( IDe Pontieu et al 



2007bl ). because of their mostly upward motion, 



short lifetimes, high upward velocities, and occurrence in a CH. 
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3. Results 

We study individual packets of transverse waves by analyzing the horizontal displace- 
ment of the medial axes of the spicules as a function of time and height (Figure [2]) . The 
middle panels show the horizontal displacements as a function of time. In both examples 
(top and bottom row), we can see oscillatory patterns as a function of time for each separate 
height (shown with color). To analyze the propagation, we focus on the peak excursions of 
the oscillations. The white circles in the middle panels indicate the location of peak excursion 
for an individual height. This peak excursion is derived separately for each height by fitting 
a Gaussian to the displacement vs. time curve for a timerange of 21 timesteps around the 
local peak in the horizontal displacement vs. time plot. We can easily see in the top- middle 
panel that several of the peak excursions move upwards (in height, i.e., color) with time 
(y-axis). This means there are upward propagating waves along this spicule. Here we define 
a set (i.e., for a range of heights) of peak excursions as one wave packet. We calculate the 
phase velocity of the wave packet at a certain height h by performing a linear fit of the peak 
excursion curve (in the time vs. height domain) over 9 heights centered around height h. @ 

The results are shown in the right panels of Figure [2] with phase speed shown in colors 
and numerals. We can clearly see in the top-right panel that one of the wave packets 
propagates upward with a phase speed that increases with height (white arrow). In the 
bottom-right panels, we see an example of the interaction between upward (at low heights) 
and downward (at greater heights) propagating wave packets. The intermediate height 
range shows evidence of a wave packet with very high phase speed. We speculate that 
this is evidence for a standing wave, as a result of the superposition of the upward and 
downward propagating packets. While a pure standing wave would have an infinite phase 
speed, we cannot expect to measure such a speed in our data, given the uncertainties in the 
measurements. In addition, if the upward and downward propagating waves do not have 
exactly the same amplitude, we would observe a partially standing wave with very high (but 
not infinite) phase speed. Here, we define waves with speeds of more than 500 km s~ x as 
standing waves, because the chromospheric Alfven speed is likely to be significantly less than 
that. The interaction of the waves is also illustrated in online animations B and C, with red 
and blue representing upward and downward propagating waves, respectively, and yellow for 
standing waves. 

Next, we study the dependence of the phase speed with respect to height and time 



2 To calculate the phase speed, we focus on individual excursions as opposed to performing a cross cor- 
relation on the whole displacement vs. time curve. The latter does not provide reliable values when a 
horizontal-time plot shows only one peak excursion. 
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(since the beginning of the spicule). Figure |3] shows the distribution of the phase speeds. 
On the horizontal axes we use the time lag between different heights, which is equal to 
(velocity) -1 , instead of the phase speed (shown on the top horizontal axis). The left panel 
indicates the height dependence of the phase speed, with the diamond symbols marking the 
median phase speed for each height range. The median for the phase speed along the whole 
length of the spicules is 270±30 km s _1 . This is roughly consistent with estimates for the 
Alfven speed in the quiet Sun, where the magnetic field strength is estimated to be 10 G 
(ITruiillo Bueno et al.l 120051 ) and the plasma number density is of order 10 10 cm -3 ( iBeckers 
1968h . The median phase speed clearly increases with height and reaches values in excess of 



1000 km s _1 at greater heights (i.e., standing waves dominate). The right panel indicates the 
dependence of the phase speed on the time elapsed since the spicule was first detected. We 
note that the start time is not the time of spicule formation, but the time at which the spicule 
fulfills our length criteria. We find a tendency for higher phase speeds to occur during the 
earlier phase, then upward waves are dominant, and finally the phase speed become higher 
again towards the end of the spicule's life. 

We also measured the wave periods and amplitudes (Figure 0J for wave packets that 
show at least two peak excursions at each height. The medians of amplitude and period are 
0.077±0.072" (55 km, 8x) and 45±30 s (P), respectively (with the standard deviation as 
error). Assuming a sinusoidal motion, we can calculate the median of the velocity amplitude 
(v = 27r5x/P) as 7.4±3.7 km s _1 . We can also estimate the velocity amplitude from a least 
squares fit of a linear function between the period and the amplitude. The best fit occurs 
for a velocity amplitude of 8.1 km s _1 . This value is consistent with the result deduced from 
the assumption of a sinusoidal oscillation. Using these parameters, and assuming a plasma 
number density of 10 10 cm -3 , we would roughly estimate the Poynting flux to be 2.5 xlO 5 
erg cm -2 s _1 if the filling factor were one. It is however difficult to assess the filling factor 
of the waves using current data. 



4. Discussion 

We investigated the lateral motions of spicules observed with Hinode/SOT, and detected 
numerous standing (20%), and upward (59%) and downward (21%) propagating waves. 

The phase speed increases significantly with height along the spicules. This could be 
caused by a combination of two factors. First, at greater heights the Alfven speed is higher 
since the plasma density is lower. If we assume that the plasma number density is roughly 
proportional to the observed intensity, the ratio of Alfven speed at every 5" compared to 
that at the surface would be 1.00, 1.99, 2.95, and 3.39 (164, 327, 484, and 556 km s" 1 with 
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164 km s" 1 at 0"-5"). Considering that the magnetic field strength is expected to decrease 
with height, the range of phase speeds at heights below 15" thus seems reasonable. However, 
the phase speeds at heights above 15" seem too high to be explained by the change in Alfven 
speed. A possible explanation is that the upward and downward propagating waves are 
more often superposed at greater heights. This would arise naturally in a scenario in which 
numerous waves are excited at lower heights, reflected at the transition region (TR) at the 
top of spicules, resulting in downward propagating waves at greater heights (just below the 
top). As a result, we would observe a combination of upward and downward propagating 
waves near the top of spicules. Given the high reflection coefficient of the TR (see, e.g., 



Hollweg et al.l Il982[ ). the downward propagating waves are expected to have amplitudes 
roughly equal to those of the upward propagating waves. However, some leakage of the wave 
energy does occur into the corona, so that the amplitudes are not quite equal. This means 
that superposition of the counterpropagating waves would lead to (partially) standing waves 
with very high phase speeds. This is exactly what we observe. We note that the spicules do 
not live very long (especially compared to the wave periods). This means that the reflected 
waves statistically will not reach back down to the bottom of the spicules. In summary, 
the combination of the increasing Alfven speed and the reflection off the top of the spicule 
provide a good explanation for the observed change in phase speed with height. It is clear 
from our observations that the observed phase speed is not necessarily equal to the local 
Alfven speed, and may in fact be very different from it because of the strong superposition 
of upward and downward propagating waves. Our results indicate that seismology of the 
spicular environment using Alfven waves requires not only very careful analysis, but may in 
fact be highly problematic and unreliable. 

This picture is confirmed when we look at the dependence of the phase speed on time: 
the phase speed decreases from the early phase to the middle phase, followed by an increase 
during the later phase of the spicule lifetime. While the start time is somewhat arbitrary, the 
definition is the same for all spicules, so it is reasonable to interpret the change as a function 
of time since spicule initiation. During the early phase, the spicules are shorter so that the TR 
is also located at lower heights. Hence, upward waves excited in the lower atmosphere would 
reach the TR sooner and are immediately reflected. This can only occur if the Alfven speed 
(~ several hundred km s _1 ) is greater than the expansion velocity of the spicule (~ 50 — 100 
km s _1 ), which seems reasonable. Therefore, most wave signals at this stage consist of both 
upward and downward propagating waves, and the observed phase speed is higher. During 
the middle phase of the spicule lifetime, upward propagating waves are dominant. This could 
arise in a dynamic situation where wave packets are continuously propagating upward from 
the bottom, with the top of the spicule moving upward, and the phase speed is determined 
over the whole length of the spicule. During the middle phase there is a lower probability of 
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the reflected waves reaching all the way back to the bottom of the spicule, so that there is 
a higher ratio of upward propagating waves. Towards the end, the reflected waves have had 
the time to propagate all the way back to the bottom (again resulting in a mix of upward 
and downward propagating waves). This scenario is speculative: the spicular environment 
is highly dynamic and the mix of waves will sensitively depend on the difference between 
the Alfven speed and the velocity of the spicule (or the TR), and on the steepness of the 
TR. Our observations should be compared with and will help const rain detailed numerical 
mode ls of Alfvenic disturbances on chromospheric spicules (see, e.g.. iMatsumoto fc Shibata 
2010h . 



Such models will also have to explain why we see that upward propagating waves are 
three times more likely than downward propagating waves. At first glance, this may b e 
seen as inconsistent with the idea that the TR has a high reflection coefficient (jDavilalll99fl ). 
Perhaps the TR is not as sharp as traditionally assumed (e.g., the hot plasma lo cated beyond 
the height of cool spicules that was recently reported by lDe Pontieu et al.ll201ll )? Or perhaps 
the dynamic spicular environment with its competition between the upward spicule velocity 
and the Alfven speed explains the dominancy of upward propagation? 

In our analysis, we have neglected the inclination of spicules. Since almost all of the 
detected spicules were within less than 10° from the vertical, the effect on the phase speed 
is small (of order 10%). We also do not know where each spicule is rooted: some are located 
in front (or behind) of the limb, while others are located right at the limb. As a result each 
height measurement is actually a mix of heights above the solar surface, so that the phase 
speed distribution at lower heights includes information from greater heights. Given that the 
distribution at greater heights includes only those heights and shows greater phase speeds 
(Figure |3]), the phase speed distribution at low heights would actually shift to the right (i.e., 
more upward propagation) if the information from greater heights were removed. In that 
case, the phase speeds would be better represented by dominance of upward propagating 
waves with phase speeds of 100-20 km s~ l . This is close to the values extrapolated to lower 



heights from the measurements by iMcIntosh et al.l (1201 ll ). 



The waves have rather short periods (< 100 s) and small amplitudes (< 0.3"). This is 
an artifact of our analysis. Because we require at least one full "wiggle" (i.e., a half period) 
during the spicule lifetime, we cannot measure the propa gation or properties of th e lower 
frequency waves (P > 100 s) that carry much more energy (IDe Pontieu et al.ll2007al ). These 
lower frequency waves are not accessible with the method used, but are present in our dataset 
as well. 



The observed high frequency waves carry a significant energy flu x of order 2.5 xlO 5 
in principle enough to play a role in heating the corona ( jWithbroe Noyes 



erg cm 2 s 1 
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19771 ). However, our observations suggest that most of these waves likely do not reach the 



corona in sufficient quantity because of the reflection off the TR. We note that low frequency 
waves are even more effic iently reflected, but carry a much larger energy flux into the corona 
(foe Pontieu et al.ll2007ah . 



Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as 
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by KAKENHI (21-8014) and Grant-in-Aid for "Institutional Program for Young Researcher 
Overseas Visits" from the Japan Society for the Promotion of Science. B.D.P. was supported 
by NASA grants NNX08AL22G and NNX08BA99G. 
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Fig. 1 — Context 17lA image (left) from SDO/AIA at 11:00 UT on 29-Jan-2011 highlighting 
the location of the Hinode/SOT field of view (top panels, in red). The right panel shows 
the sum of 10 consecutive AIA images taken during 2 minutes which has been treated with 
unsharp masking to bring out the faint structures off-limb. Notice the large helmet streamer 
like configuration which has its roots to the north of the SOT field of view. This indicates 
that SOT was most likely pointed at a region that is part of the southern polar CH (i.e., 
open field conditions). The bottom panel shows a simultaneous Hinode/SOT Ca II H image 
illustrating the prevalence of chromospheric spicules. 
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time=83.2 sec 





Fig. 2.— Two examples of detected spicules and their transverse displacement. Left: 
Snapshot. The detected spicule is highlighted with a colored line (color representing height). 
Center: Time variation of horizontal displacements at each height of the spicule. The color 
represents height (same scale as left panel). Peak excursions of each wiggle are marked with 
white circular symbols. Right: Peak excursions and the phase speeds in a time-height map. 
The numerals denote the values of the phase speed (km s _1 ). The upper panel shows that 
most of the waves propagate upward. One of the wave packets is seen to propagate with a 
speed that increases with height. The lower panel shows that there are a variety of wave 
packets along the spicule. In one case a wave packet propagates upwards from low heights 
and encounters a downward propagating wave packet (greater heights). A standing wave is 
produced in the middle as a result of superposition of the propagating waves. 
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Fig. 3. — Height and time dependence of the measured phase speed distributions. Left: 
Histogram of phase speed binned for different height ranges (separated by five arcseconds) 
above the solar surface. The diamonds indicate the median values of the phase speed at 
each height range. Right: Same as the left panel, but now for different time periods since 
the first detection of the spicules (binned every 40 seconds). 
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Fig. 4. — Histogram of the measured amplitudes and periods (top row), and velocity am- 
plitude (bottom left). Dashed lines show the median value. The velocity amplitude is deter- 
mined from the amplitude and period by asssuming sinusoidal oscillations. The bottom-right 
panel shows a scatter plot of period and amplitude. The full line is a best fit which shows 
that a velocity amplitude of about 8 km s _1 fits best. This value is consistent with the me- 
dian of the velocity amplitude deduced from the sinusoidal oscillation assumption (bottom 
left). 



